Theory of Spin Susceptibility in Frustrated Layered Antiferromagnets 
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The self-consistent treatment of real and imaginary renormalizations in the dynamic spin suscep- 
tibility x{c{,Lj) for the frustrated Heisenberg model reproduces for cuprates at low doping: a spin 
spectrum ojq, a saddle point for q « (7r/2,7r/2), nearly constant g-integrated susceptibility X2d(<^) 
for u < 150 meV and a scaling law for X2d{i-o)- Frustration increase (optimally doped case) leads 
to a stripe scenario with an Li^q-saddle point at q ss(7r; n/2) and X2d(i^) peak at w « SOmeV. The 
obtained x(q, <^) describes neutron scattering results and leads to well-known temperature transport 
anomalies in doped cuprates. 

PACS numbers: 74.72.-h, 71.27.+a, 75.20.-g 
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The investigation of the dynamic spin susceptibility 
x(q, uj) is a key problem for understanding the physics 
of layered high- Tc superconductors (HTSC) in both low 
and optimally doped regimes. The inelastic neutron scat- 
tering (INS) measurements in cuprates 0-0 re- 
vealed a sharp resonant magnetic excitations peak of 
X (q, t^) = I^xi^i^) which corresponds to the anifer- 
romagnetic (AFM) wave vector Q = (tt, tt) at a reso- 
nant energy Er ~ 30 meV and a low-temperature peak 
at close energies for g-integrated susceptibility X2Di(^, T), 
X2d{'^, T) — f dqx (Qi ^1 T). At low doping INS demon- 
strates the scaling of magnetic response ~ the universal 
law for X2d(w,T) 4] which states 



X2d(^,T) 

X2d{^,T 0) 



(1) 



In the mentioned regime the spin excitation dispersion 
a;(q) was measured across the entire Brillouin zone ^3 
and it was found that aj(q) is anisotropic around the 
magnetic zone boundary (a saddle point at q ~(7r, 7r/2)). 

The aim of the work is to present a theory for the dy- 
namic spin susceptibility x(q, a;) within the frustrated 
S = 1/2 Heisenberg model taking into account real 
and imaginary renormalizations extracted from the ir- 
reducible Green's function M(q,(jj) so as to describe the 
mentioned experimental results in both doping regimes 
in the framework of one self-consistent approach. Our 
analysis is based on a spherically-symmetric treatment 
of the spin system which was introduced by Shimahara 
and Takada Q and generalized in Q. 

The recent microscopic theoretical progress in the in- 
vestigations of x(q, w) 0, IE is based on t — J model 
treated within the memory function approach (MFA). 
This approach demonstrates x (Q: i^) peaks and the scal- 
ing law. It is close to our treatment but it has difficulties 
in an analytical calculation of the explicit form for a;(q) 
and a spin gap and as a result in the self-consistency 
procedure. Relative to MFA our theory gives such new 



results as a demonstration of w(q) - saddle point at 
q w(7r;7r/2) and a new analytical form for a scaling law 
for small frustration case (strongly underdoped regime). 
For the large frustration (the regime close to optimal dop- 
ing) we reproduce not only x '(Q,cj) peaks, but also the 
peaks of X2d('^) demonstrating a stripe scenario (see 
for a review). For the latter case we calculate also the 
resistivity p{T) and the Hall coefficient R{T) in order to 
be sure that the found x(Qi reproduces the well-known 
temperature anomalies in kinetics. 

The Hamiltonian of the model has the form 



— > — > 1 X ^ — > — !■ 

s i s i+g + -I2 s i s i+d 



(2) 



It describes the frustrated system of localized 5=1/2 
spins on a square lattice, where Ii is an AFM interaction 
constant for nearest, I2 ~ for next-nearest neighbors, g, d 

- vectors of nearest and next-nearest neighbors. We use 
standard variable p ("frustration parameter") p = I2I I, 
Ii = (1 — p)I, I2 = pi, I = I1+ I2 as a measure of frus- 
tration, hereinafter we treat all the energetic parameters 
in the units of / and put / = 1. We suppose that the 
frustration (term I2) simulates the influence of doping. 

Following IE la ^ we calculate x(q;^) — 

— ((S'q I S'"q))aj - two-time retarded Green's function 
by the irreducible Green's function method. The dy- 
namic spin susceptibility can be written as x{^t^) = 
—Fq/{uj'^ — ojq — M(q, cj)), where ojq is the spin excita- 
tions spectrum, M(q, ui) = M +iM" - Fourier-transform 
of a new complex three-site irreducible retarded Green's 
function, its analytical properties are the same as those 

of x(q,w). 

Spectrum ujq and the numerator Fq have a cumber- 
some form but they are expressed explicitly over five 

g+d 

This 
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spin-spin correlation functions Cg, Cd, Cg+g, Cd+d, c 
i, Cr - {S^S^+,) = (2^)-2/dqcqe^'i' 
allows to write down and solve numerically sclf 
consistent system through the usual relations c, 
(S^Si^) = -TT-^ JdLjnB{u;)Im{{S^ \ S%))^+,S: 
]\[-i/2 e^^'^'^S^. The set of equations includes also the 
q 

sum-rule condition Cr=o = 1/4. The system is solved at 
every fixed T and p. 



q ~ 

ca 
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The imaginary part M"(qL w) is an odd function of 
u>. In the simplest approach [g] we put M"(q, t^;) = — a;7, 
where the damping 7 is taken to be independent on q and 
Lo. We take the real part M as M ~ |sin(ga;) sin((7j^)|'^ 
and introduce a renormalized spectrum tjq = Wq + 

(A| sin^a; sinQj^l)'^. The choice of M functional form is 
dictated by the condition that M represents the square 
harmonic different from those involved in the functional 
form of ti'^(q). Though the A-renormalization is zero 
along the lines T - X and X - M (F =(0, 0), X =(0, tt), 
M =(7r, tt)) and mainly modifies the top of the spectrum, 
it changes the spin gap Am = ujq due to self-consistency 
of calculations. So the dynamic spin susceptibility 

— F 

X(q,w)==^ — , (3) 

contains two parameters 7 and A. 

We relate the dielectric limit to the case of extremely 
small frustration p = 0.04. In the inset of Fig.l the 
spectrum S3(q) is presented in this limit for T = 0.1, 
7 = 0.025 and A = -1.0 (T - lOOA' for / - 100 mel^). 
The spectrum is almost linear on q — |q — Q| up to 
LUo ^ 1-5. It can be found that for fixed q there is a 
well-defined x(q, uj) peak on uj which is related to the 
spectrum cD(q). More exactly, the maximum of x(q, w) 
on UJ corresponds to the frequency close to w(q), but al- 
ways a bit smaller (due to damping 7). For / = 1.2 meV 
a spin- wave velocity he « 900 meV A is close to the 
value given in j^H- As it is seen from the inset of Fig.l, 
in accordance with the experiments 13J, the dispersion 
tD(q) is anisotropic around the magnetic zone boundary 
and has a saddle point close to q = Q/2 (u;(q =(0, tt) > 
a;(q =(7r/2, 7r/2))). Note that in contrast to our treat- 
ment one needs to adopt a ferromagnetic second-neighbor 
exchange I2 with p < —0.1 for the explanation of such 
an anisotropy in the framework of the linear spin-wave 
theory [3- 

In Fig.l X2D(t^) is given for p = 0.04, 7 — 0.25T in 
two cases: T = 0.1, A = -1.0 and T = 0.3, A = 2.0^/^. 
The A-values are chosen from the condition that the re- 
sulting spin gap should be approximately linear on T : 
Am(T = 0.1) = 0.048, Am{T = 0.3) = 0.134. Below we 
show analytically that in the low-frustration limit this 
linearity is the necessary condition for the scaling law. 
The qualitative coincidence of calculated function with 
the experiment 0, is seen - X2d(w) is nearly con- 
stant in a large oj interval and increases for oj > 150 meV. 

Now we treat the scaling condition which leads to 
a strong limitation of j(T) dependence. Fig. 2 rep- 
resents the scaling functions f{uj/T). Solid lines b 
and c correspond to temperatures T = 0.1 and T — 
0.3 respectively and are calculated for the same pa- 
rameters, as in Fig.l (that is, in particular, for 7 — 
0.25T). The dash-dotted line a is the best fit for ex- 
perimental scaling in Lai,geSro.o4Cu04 Q fex{ijj/T) = 
(2/7r)arctan{0.43(w/T) -I- 10.5{uj/T)^}. It is approxi- 
mately a step function on {uj/T) smeared through S = 




FIG. 1: X2d{uj) for frustration p = 0.04: T/I = 0.1 - solid, 
T/I — 0.3 - dashed curve (damping 7 — 0.5T). Inset: self- 
consistent spectrum c5(q) for p — 0.04 and T/I = 0.1 
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FIG. 2: Scaling curves /(cj/T) for p = 0.04: the dashed- 
dotted line a - best fit for scaling in Lai.g6Sro.o4Cu04 Q; 
solid lines b and c were calculated at T/I = 0.1; 0.3 for 
7 = 0.25T; dotted line e - scaling law Q - see text; solid 
thin curve d- destroyed scaling T = 0.3 and 7// — 0.3 > 
0.25T/I = 0.075. Inset: Scaling for p = 0.1: solid curves 
f,g,h respectively for T/I = 0.1;0.2;0.4, 7 = T. The solid 
thin curve 1 - T/I = 0.4 and 7// = 0.1 < T/I = 0.4. 



A(cj/T) ~ 0.25. The calculated curves have close value of 
6. Note that the value of S strongly restricts the 7(T) de- 
pendence. For example, thin curve d in Fig. 2 corresponds 
to T = 0.3, A = 2.01/3 and 7 = 0.3 > 0.25r = 0.075. As 
a result curve d strongly deviates from curves a ( fex ) 
and c (/T=o.3,7=o.25r) and it has S S> 0.25. 

The analogous picture for the frustration p = 0.1 is 
shown in the inset of Fig. 2 for 7 = T and A = We 
relate this case to strongly underdoped F-cuprates [l^ . 
The calculated scaling functions are given for the tem- 
peratures T = 0.1,0.2,0.4 by the solid curves f, g, h 
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respectively. It can be seen that they are close to the 
experimental /^.^{uj/T) — (;|) arctan{a(a;/T)}, a ~ 1 
hj. The thin curve 1 corresponds to T = 0.4 and 
7 = 0.1 < T = 0.4. Comparison of the curves h and 
1 explicitly demonstrates that, as in the previous case 
p = 0.04 , the scaling is destroyed when 7(r) deviates 
from a linear law. 

Thus, the above results demonstrate that the scaling 
law holds if 7 is a linear function on T. In the limit 
p <C 1 this point can be clarified analytically taking 
/mx(q, uj) = 7cji\j/{(tj^— a;q)^+7^a;^} It is obvious from 
the inset of Fig.l that for uj/T < 2, T/I < 0.3 the main 
input to X2d(w) is given by the region q < q^; cqo ^ I is 
the largest energetic parameter. Thenw^(q) R; 
and simple integration gives for uj < cqo 



A2 - 



c q 



X2d{uj) 



47rc2 



$(tj,AM,7); fore 
7r + $(w, Am,7); for- 



< 1 
' > 1 



(4) 



arctan 



9^ 



72^2 + (c2q^ + A2^ 



Cj2) 



Here Fq is the averaged smooth function Fq. 

In the limit under consideration the scaling denomina- 
tor X2d{^,T — !■ 0) is almost constant in a wide w-range 
and scaling is ruled by X2d(w,T). Accepting in lin- 
ear 7 = aT and Am = I3T one obviously gets the scaling 
{X2d{'^,T) becomes X2d{^/T)). So in the mentioned 
approximations the scaling function can be written as 



-f arctan 



a{uj/T) 



(/32 - (cu/m 



(5) 



In contrast to numerous experimental fittings by simple 
arctan, the scaling function f{uj/T) Q is described by 
'switched' arctan law and contains a microscopic infor- 
mation on Am and 7. The switching by 0-function takes 
place at w = Am- 

In Fig. 2 / is represented for a — 0.25, /3 = 0.5 by 
dotted line e and it coincides with /t=o. 2, 7=0.05- The 
scaling function /(w/T) with slightly different parame- 
ters a = 0.25, /9 = 0.43 is almost indistinguishable from 
experimental fex{'-^/T) Let us remind, that in the 
above calculations we have taken 7 ~ T and such A(T) 
that the self-consistent calculations led to A ~ T. 

So in the dielectric limit (small p) the model leads to an 
adequate description of experimental results. The scaling 
law strongly restricts 7(T') dependence. 

Now we turn to the case p = 0.28 which corresponds to 
Am ~ Ax = '^q=x- We relate this case to the optimal 
doping. Calculated data are presented for T — 0.025 and 
0.05 with 7 = 0.38-I-0.8T (in contrast to low frustration 
limit here 7 does not tend to zero at T 0) and A = 
10.0^/3. For T = 0.025 and 0.05 the gaps are Am = 
0.197, Ax = 0.179 and Am = 0.228, Ax = 0.210. 

Fig. 3 shows the Q-peaks, i.e. x (Qi'^)^ for T = 0.025 
and T = 0.05. They are also in good agreement with 
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FIG. 3: X (Q,^) for T/I = 0.025 (solid) and T/I = 
0.05 (dashed curve) Inset: x^oiuj) - solid line and 
X2D(tt')(2n_gose-l-l) - dashed line for T /I — 0.05; Lower inset: 
The spectrum u^/I for p = 0.28, T/I = 0.05. 

the experiment In the inset of Fig. 3 the calculated 
spectrum Wq is shown for T = 0.05. The dispersion 
w(q) has the following new features: the saddle points 
close to q =(7r; 7r/2); (7r/2; tt) and ui{q) changes weakly 
along X — M direction. That is why X2d('-^) has a peak 
at a; > 2 3 Am- This is explicitly seen in the an- 
other inset of Fig. 3 which gives X2d(w) (solid line) and 
X2D(w)(2nBose + 1) (dashed line) for T = 0.05. These 
curves qualitatively correspond to the experimental ones 
[m for optimally doped curates. It is clear that the 
shown behavior of w(q) and X2d{^) is a result of a stripe 
scenario if we remind that the increase oip drives the sys- 
tem to a state which is close to a coherent superposition 
of two semiclassical stripe phases with Ax = Q - It can 
be shown that Cq — (S^St^) is qualitatively different for 
small and large frustrations- For p < 0.1 the structure 
factor Cq has an extremely narrow peak at q = Q. For 
p — 0.28 the structure factor has peaks at q = Q = M 
and at q = X. With p increase the peaks at X points 
increase and the M -peak disappears. 

We capture this physics taking a spin-only model. 
But this model is too simple to reflect a well-known 
low-energy incommensurate magnetic excitations at wave 
vectors close to Q at optimal doping. It is obvious that 
one needs to introduce explicitly the spin-hole scattering 
to describe this feature. 

To check the applicability of the obtained spin suscep- 
tibility x(q, 1^) for the kinetics of the optimally doped 
HTSC we calculate the in-plane resistivity p{T) and the 
Hall coefficient Rh{T) in the framework of the spin- 
fermion model with the Hamiltonian 

H ^ Ho + Hi (6) 

k,<T ^ k,q,7i,72 

The hole spectrum £k is obtained from the calculation of 
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FIG. 4: The resistivity p{T/I) and cotangent of the Hall angle 
coteniT^/I^) at 10 T (upper mset). The solid lines - for 
x(q, t^) ©• The dashed lines - for overdamped Xojjd(q, <^)- 
Lower inset: the spectrum ek (in eV) given by the curves 
£k = const; bold curve - Fermi line for optimal doping. 



the lower spin-polaron band in a six pole approximation 
[l7| and is shown in the inset of Fig. 4. 

It is well known that scattering by the spin- 
fluctuations with momenta Q leads to a strongly 
T-dependent anisotropy. To take it into account 
the equation of motion for the non-equilibrium den- 
sity matrix p'^^^ = exp(— iJo/T)-F is solved 
by seven-moment approach F — '^i^i^rVi^h 

= Z]k,o--^'(k)aL«ka- The moments F/(k) are 
taken to be polynomials in velocity components 
Vk — dek/hdk and their derivatives: Ff'Ck) ~ 

U,x („V\2.,x dvg y dvl X dv^^„,x ( x\3 8<„,xl 
l^k' V^J ^^k' dy ^k' dy ^k' dx dy ^k' l^kJ ' dx ^k/- 



The detailed expressions for p{T) and Rh{T) are given 
in jl^. The susceptibility x{^.^^) ^ is involved in 
scattering integrals. To clarify the importance of the 
form © we also present the results for widely used 
so-called overdamped susceptibility Yr>„rf(q, a;) (when uj^ 
term in the denominator of JSJ is omitted) [lil uX\ . 

The results presented in Fig. 4 are obtained for p — 
0.28, / = 100 meV and J = 200 meV. The plots 
are the resistivity p{T) and the Hall angle cotangent 
cotQu — Pxxl {RhB) (in the inset) obtained for the 
x(q, Lo) ^ - solid lines and for Yo„rf(q, ui) - dashed hues. 
In accordance with the experiment |2l|. the p{T) curve 
exhibits a temperature dependence close to a linear one 
starting from low T with the value p{AOQK) / p{\QQK) 
~ 5. It can be shown that Xovd{f\,^) approximation un- 
derestimates the scattering for large As a result at 
p(T)ovd < p{T) and, as it is seen from Fig 4, in some 
temperature regions p{T)ovd has a different curvature. 
The cot Qh exhibits nearly linear behavior on in a 
wide temperature range, however, at low temperatures 
deviation from linearity is obvious. It seems hopeful that 
the self-consistent spin susceptibility x(q, ^) allows to de- 
scribe experimental temperature anomalies of two kinetic 
coefficients simultaneously. 

In summary, we have made a systematic self-consistent 
study of the spin problem in 2D frustrated Heisenberg an- 
tiferromagnet. Key features of the model - temperature 
dependence of the damping in low frustration limit and 
the appearance of saddle points of the dispersion lo{c\} 
close to q =(7r;7r/2); (7r/2;7r) in the case of strong frus- 
tration increase - allow to relate the results to a wide 
hole doping interval in cuprates. 

Work was supported by Russian Fund of Fundamental 
Investigations and Russian Science Support Foundation. 
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